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We study the properties of pions in twisted mass lattice QCD (with two degenerate flavors) using 
chiral perturbation theory (%PT). We work to next-to-leading order (NLO) in a power counting 
scheme in which m q ~ aAq CD , with m, the physical quark mass and a the lattice spacing. We 
argue that automatic O(a) improvement of physical quantities at maximal twist, which has been 
demonstrated in general if m q 2> aAq CD , holds even if m q ~ aAq CD , as long as one uses an 
appropriate non-perturbative definition of the twist angle, with the caveat that we have shown 
this only through NLO in our chiral expansion. We demonstrate this with explicit calculations, for 
' arbitrary twist angle, of all pionic quantities that involve no more than a single pion in the initial and 

t final states: masses, decay constants, form factors and condensates, as well as the differences between 

alternate definitions of twist angle. We also calculate the axial and pseudoscalar form factors of the 
pion, quantities which violate flavor and parity, and which vanish in the continuum limit. These 
are of interest because they are not automatically O(a) improved at maximal twist. They allow a 
determination of the unknown low energy constants introduced by discretization errors, and provide 
ly-^ tests of the accuracy of xPT at NLO. We extend our results into the regime where m q ~ a 2 Aq CD , 

and argue in favor of a recent proposal that automatic O(a) improvement at maximal twist remains 
valid in this regime. 
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INTRODUCTION 

Twisted mass lattice QCD (tmLQCD) is an alternative regularization for lattice QCD that has recently 
received considerable attention. 1 It has the potential to match the attractive features of improved staggered fermions 
(efficient simulations absence of "exceptional configurations" 0, 0(a) improvement at maximal twist 5], operator 
mixing as in the continuum while not sharing the disadvantage of needing to take roots of the determinant to 

remove unwanted degrees of freedom. Initial numerical investigations of the unquenched theory 0,0 have, however, 
found surprisingly large discretization errors, which manifest themselves as a non-trivial phase structure. 2 While these 
errors do not pose a fundamental problem, and, indeed, appear to conform to the expectations of chiral perturbation 
theory applied to the lattice theory |12l Ha . [lj, , they do suggest that a thorough investigation of the impact of 
' discretization errors on the theory is called for. 

Such an investigation is possible by applying the methods of chiral perturbation theory to tmLQCD at non-zero 
lattice spacing. The discretization errors are included systematically in a joint expansion in the lattice spacing, a, 
and the quark mass, m q . The resulting "twisted mass chiral perturbation theory" (tmyPT) has been formulated 
previously 0[H|ig, 

building on earlier work for the untwisted Wilson theory [T3 . [T^ . fl8| . It has been used to study 
pion masses and decay constants in the regime where m q ^> o,Kjq CD [l^, and to study the phase structure of tmLQCD 
when m q ~ a 2 h^ CD 0, El El El El ■ Here we build upon our previous work [l^, in which we determined the 
chiral Lagrangian at next-to-leading order (NLO) in a power counting scheme in which we treat m q ~ clAq CD . We use 
this to study all the quantities involving pions that do not involve final state interactions and are thus straightforward 
to calculable in simulations: masses (previously calculated in this regime in Ref. |14|). condensates, vacuum to pion 
matrix elements, and matrix elements between single pion states (i.e. form factors). The operators we use are the 
vector and axial currents and the scalar and pseudoscalar densities. 
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1 For a recent review sec Rcf . 3 1 . 

2 Indications of surprisingly large discretization errors have also been fou nd in quenched tmLQCD at maximal twist llll . On the other 
hand, use of an improved gauge action reduces the discretization errors ITU . 
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Calculating these quantities allows us to address some of the issues that are of concern when using tmLQCD. In 
particular, what is the impact of the breaking of flavor and parity symmetries at non-zero lattice spacing? What 
are good quantities to use to determine the typical size of discretization errors? How does the automatic 0(a) 
improvement at maximal twist manifest itself in tm%PT? How do different definitions of the twist angle, all of which 
are equivalent in the continuum limit, differ at non-zero lattice spacing? How can we test the reliability of tm^PT 
applied at a given order? We return to these questions in the concluding section. 

Another issue, which is of particular importance for practical simulations, is the condition on the quark mass that 
must be enforced in order that physical quantities are automatically 0(a) improved at maximal twist. Does one need 
to enforce (A) m q ^> aAg CD , (B) m q ~ aAq CD 3> a 2 Aq CD , or (C) m q > ca 2 Aq CD , with c a determinable constant of 
0(1)? If one takes or 1 = 2 GeV and A QC d = 0.3 GeV, these conditions become (A) m q > 45 MeV, (B) m q > 7 MeV, 
and (C) m q > 7 MeV, respectively. Since the average light quark mass is ~ 3 MeV (at a renormalization scale of 
2 GeV — for a recent review see Ref. H^), we need to be able to use condition (B) if lattice masses are to allow 
extrapolation to physical values, or condition (C) if they are to approach them. Condition (A) is very restrictive. 

On the theoretical side, Refs. 0,0 argue that, in general, one must use the strongest condition, (A), although it 
is possible that the intermediate condition, (B), can be used if one is dealing with O(a) improved quantities, e.g. if 
using an 0(a) improved quark action. By contrast, it has been recently proposed in Ref. |2l| that there is no lower 
limit on m q in order for automatic 0(a) improvement to apply, as long as one uses an appropriate definition of twist 
angle. 

We will argue here in favor of the weakest condition, (C), with the constant c either of O(l), or vanishing, depending 
on the sign of an unknown constant in the chiral Lagrangian. Our arguments are based on tnr^PT applied in the 
regimes where conditions (B) and (C) hold. The former we call the "generic small mass" (GSM) regime, the latter 
the "Aoki" regime (since it is in this region that non-trivial phase structure due to discretization effects appears |23j). 
For most of the paper we consider the GSM regime, in which we expect the bulk of future simulations to be done. We 
give a general argument for automatic 0(a) improvement in this regime (which requires the use of an appropriate non- 
perturbative definition of twist angle), and support it with results for the physical quantities listed above. We then 
extend the results to the Aoki regime, where, working to leading order (LO) in an expansion in which m q ~ a 2 AQ CD , 
we argue that condition (B) can be relaxed to (C). 

We stress that our discussion of O(a) improvement is within the context of chiral perturbation theory up to a given 
order (NLO in the GSM regime, LO in the Aoki regime). This means that we only control those 0(a) corrections that 
are expected to be dominant in the regimes under study. For example, in both regimes we control corrections to the 
pion decay constant, f v , which are of relative size oAqcd, but do not control those of relative size am q . The latter 
corrections arise from terms in the chiral Lagrangian of 0(a 3 ) and 0(a 4 ) in the GSM and Aoki regimes, respectively, 
which are, in both cases, of higher order than we work. 

The remainder of this article is organized as follows. In the next section, we recall the definition of tmLQCD and 
the construction of the corresponding NLO chiral Lagrangian, and use this to discuss the condition on the quark mass 
needed to obtain automatic 0(a) improvement at maximal twist. In Sec. lllll we use the Lagrangian, including sources 
for currents and densities, to determine the twist angle non-perturbatively, and to calculate the physical quantities 
listed above. In the final subsection, we give a summary of our results and a suggestion for how to use them in 
practice. In Sec. IIVI we show how our results extend into the Aoki regime, and comment on the considerations of 
Ref. 2jJ. We conclude in Sec.E] and suggest various avenues for future work. 

Some preliminary results from this paper have been presented in Ref. 15]. 



II. THE EFFECTIVE CHIRAL LAGRANGIAN 



The theory we consider is tmLQCD with a degenerate doublet of quarks. The fermionic part of the Euclidean 
lattice action of the theory in the so-called "twisted basis" is |j, La : 

Sp = J2 ft*) [\ >( V M + V ^ - + m ° + MX), (1) 

X /I /I 

where ipi and -ipi are the dimensionless bare lattice fields (with "I" standing for lattice and not indicating left-handed), 
and V M and V* are the usual covariant forward and backward dimensionless lattice derivatives, respectively. The 
field ipi is a flavor doublet, and T3 acts in flavor space and is normalized so that rf = 1. The bare normal mass, mo, 
and the bare twisted mass, fiQ, both of which are dimensionless, are taken to be proportional to the identity matrix 
in flavor space. 

The effective continuum chiral theory is derived using the two-step procedure of Ref. [l2|. In Ref. 15], we have 
carried this procedure out to NLO in an expansion in which we treat quark masses and the leading discretization 
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errors as symmetry breaking parameters of the same size. Here we recall only the essential details. 

Following the program of Symanzik p4| , we first write down an effective continuum Lagrangian at the quark level 
which describes the long distance physics of the underlying lattice theory. Its form is constrained by the symmetries 
of the lattice theory to be 

£ ff = C g + i>{Jp+ m + i75T 3 /i)?/> + biatpia^F^ip + 0(a 2 ) . (2) 
Here C g is the continuum gluon Lagrangian, m is the physical quark mass, defined in the usual way by 

m = Z m (m - fh c )/a, (3) 

and n is the physical twisted mass 

A* = Z^no/a = Z p 1 fi () /a , (4) 

with Zp the matching factor for the pseudoscalar density. The quantity m c is the critical mass, aside from an 0(a) 
shift. This shift, and methods for determining the critical mass after its inclusion, are discussed below. Note that the 
lattice symmetries forbid additive renormalization of [io [2j. Although the 0(a 2 ) terms in Eq. J2J) are of an appropriate 
size to be included at NLO, they do not break the continuum symmetries any further than the terms explicitly shown, 
and thus do not lead to any additional operators in the effective chiral theory |l5j| . Thus we do not need their explicit 
form. 

We reiterate here that we have dropped from C e s terms proportional to a/j^ipip, arr^ipip an d a^mipi^T3ip. These 
terms are allowed by the lattice symmetries |25|. and do lead to corrections linear in a, but are of next- to-next-to- 
leading order (NNLO) in our power counting. 

In the continuum limit [where the b\ and higher order terms in @ vanish] the apparent flavor-parity breaking due 
to [i is misleading, since it can be rotated away by a non-anomalous axial rotation. Thus tmLQCD is, in this limit, 
equivalent to QCD with two degenerate quarks of mass \J m 2 + /i 2 . This has been established in detail in Ref. 0. 
In this paper, we work away from the continuum limit, taking m ~ /i ~ p 2 ~ gAq CD . Thus the symmetry breaking 
induced by the b\ term will play a crucial role. 

The next step is to match the continuum effective Lagrangian @ onto a generalized chiral Lagrangian. Working 
at NLO in our power counting, we found |15| (after simplifying using the properties of SU(2) matrices) 

C x = Ltt(D^D^) - ^Tr( X tvj + tf x ) - Ltt(A^ + £U) 

- L^iD^D^f - L 2 Tt(D^D^)Tt(D^D u ^) 
+ L 45 Tr(^S t ^S)Tr(x t S + S+ X ) 

+ L 5 {Tr [(D^D^)( X ^ + £ f x)] - Tr(D^ D^)Tr( X ^ + Y? X )/2} 
-L 68 [Tr( X ts + S t x) ] 2 _ i8 | Tr [ (x t s + s t x) 2] _ [Tr( x t E + E t x) ] 2 /2 } 

+ iL g Ti(L^D^D^ + R pu D^D v T,) + L w Ti{L^R^) 

+ W i5 Tr(D^D^)Tr(A^ + E f A) - W 68 Tr(x f Z + E^Tr^E + E f i) 

- [Tr(itS + E*i)] 2 + W w Tt(D^D^ + D^D^A) 

+ H^L^L^ + R^R^) - H 2 Trtfx) - H 2 Tr(A^ X + X f i) - H 3 Tr(A t A) . (5) 
As usual, the field £ is SU(2) matrix-valued, and transforms under the chiral group SU(2)l x SU(2)r as 

E — > LT,R\ L e SU(2) l , Re SU{2) r . (6) 

The quantities x an d A are spurions for the quark masses and discretization errors, respectively. At the end of the 
analysis they are set to the constant values 

X — > 2Bo(m + it 3 {i) = m + iT 3 ji, A — > 2W a = a , (7) 

where Bq and Wq are unknown dimensionful constants, and we have defined the quantities m, ft and a. We have 
included sources for left- and right-handed currents in the standard way using covariant derivatives and associated 
field strengths, e.g. 

£> M E = 0„E - i^E + iEr M , L„ v = d^ v - dj p + i[£„, l v \ . (8) 
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Sources for the scalar and pseudoscalar densities are similarly included by writing x = 2i?o(s + ip). Note that t^, r^, 
s and p are all hermitian matrix fields. Finally, we note that at NLO we should also include the Wess-Zumino-Witten 
term |2f|. We do not, however, consider here any processes to which it contributes (e.g. it — > 77), and so do not 
write it out explicitly. 

The chiral Lagrangian contains a number of parameters which are not determined by symmetries — the low energy 
constants (LECs). At leading order there are the usual / (normalized so that f v — 93 MeV) and Bq, supplemented 
by the additional constant Wo introduced by discretization errors [l^]- At NLO there are significantly more LECs. 
The Li are the usual Gasser-Leutwyler constants of continuum xPT (2?]]. We have introduced the combinations 
£45 = L 4 + L5/2 and = L 6 + L s 12 which are useful for two flavors. The W% are introduced by discretization errors 
at linear order in the lattice spacing 0, while the W[ appear at quadratic order 0. Again, certain combinations 
are useful: W i5 = W± + W 5 /2, W 68 = W 6 + W 8 /2 and W$ s = Wq + Wg/2. Finally, the H l involve sources alone, and 
give rise to contact terms in correlation functions of currents and densities, as well as contributing to their vacuum 
expectation values. 3 The H[ are the extra such terms introduced by discretization errors (which we did not include 
explicitly in Ref. but will need here). 

It is useful to understand the relation of our result for C x to that in Ref. In that work, the authors used 

a different power-counting, namely oAqcd ~ {m 2 + / i2 )/Aq CD . They worked to linear order in a, and without 
sources. Our results agree with theirs to the order they worked. They also changed variables using an axial rotation 
so that x was diagonal. This has the advantage of moving the twist entirely into the term caused by discretization 
errors, A — ► a(m — i^Ts)/ y m 2 + /1 2 . It was then straightforward to see that discretization errors in and are 
proportional to am/ \/m 2 + /1 2 and thus vanish under mass averaging (m «-> — m) and at maximal twist (m = 0). 

In our power counting, where discretization errors are superficially as large as quark mass effects, we must follow 
a different strategy to see automatic 0(a) improvement. The key point, first noted in Ref. |12j . is that the leading 
discretization effect [the third term in C x in JSJ] has the same form as the leading mass term [the second term in C x ]. 
This follows directly from the fact that the b\ term in the quark-level effective Lagrangian J2J) has the same chiral 
transformation properties as the quark mass term. Thus the LO chiral Lagrangian is unchanged from its continuum 
form if one uses the shifted spurion 

X' = X + A. (9) 

This corresponds at the quark level to a redefinition of the untwisted component of the quark mass from m to 

rri = m + aWo/B = (m + a)/(2B ) . (10) 

This shift corresponds to an 0(a) correction to the critical mass, so that it becomes 

m c =m c - aW Q /B . (II) 

In practice one automatically determines m c (and thus automatically uses m! rather than m) if one uses an appropriate 
non-perturbative definition of twist angle. This is discussed in detail in the following section and summarized in the 
last subsection of Sec. IIIII 

Since the LO Lagrangian takes the continuum form, the vacuum expectation value of E at this order is that which 
cancels out the twist in the shifted mass matrix: 

(0|E|0) LO = E = ™ + ° + *Ar 3 ^ cxp( ^ 0T3) ) (12) 



M 1 



where 



M' = \ X '\ = y/(m + a) 2 + A 2 (13) 
is the LO result for the pion mass-squared. If we define the physical quark mass by 



m q = \Jm! 2 + y? , (14) 

then it follows from l|12|) that 

cq = coscjQ = m'/rriq , so = sinwo = M/ m g ■ (15) 



3 Strictly speaking the Hi and H[ are not low-energy constants, since they absorb short-distance divergences. In fact, in Ref. l27l they 
are denoted "high-energy constants". Nevertheless, for brevity, we refer to them also as LECs. 
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Equation i112[l is one of the definitions of twist angle that we will use, although it is not the most simple to determine 
through simulations. It is important to realize that it differs by O(l) from the naive definition tanw = /Lt/m when 
to ~ aAq CD . When working in this mass regime, it is thus crucial to use the shifted mass to', or, equivalently at 
NLO, one of the definitions of twist angle discussed below. This point has also been emphasized in Ref. [2l| . 

We can now present our argument for why automatic 0(a) improvement holds even when m ~ /i ~ aAq CD . This 
will be borne out later by our detailed results. We begin by expressing the chiral Lagrangian in terms of x' instead 
of x- 

C x = ^Tr(D^D^) - ^Tr( X 'tE + E^) 

- LiTr^EZ^St) 2 - L 2 Tr( J D /i E J D v Et)Tr( J D^E£> 1/ Et) 

+ L 45 Tr(£> M Et£) M E)Tr(x' t E + tf X ') ~ L 68 [Tr^E + SV)] 2 

+ L 5 {Tr [(£> M Et£) M E)( X tE + EV)] - Tr(^St^E)Tr( X ' t E + £V)/2} 

- L 7 [Tr( X 'tE - Et x ')] 2 - L 8 {Tr^E + EV) 2 ] - [Tr( X 't E+E t x ')] 2 /2 } 
+ tL g Ti{L^D^D^ + R^D^DuZ) 

+ W Tr(Z) At E t Z) At E)Tr(i t E + £ f i) - W Tr(x' f E + E t x')Tr(i t E + E f i) 

- W [Tr(its + EU)] 2 + W 10 Tr^i^E + D^D^A) 

- ff 2 TSr(xV) - #'Tr(iV + X n A) . (16) 

We have dropped the terms proportional to L\q, Hi and H3, since these lead only to contact terms in correlation 
functions, which we do not need below. We have also introduced the useful combinations 

W = W 45 - L 45 , W = W 68 ~ 2L 68 , W 1 = - W 68 + L 68 , H' = H' 2 -H 2 . (17) 

This Lagrangian has the same form as the NLO continuum chiral Lagrangian aside from the "W" and H' terms. 

We are interested in the terms which are linear in a. Setting aside the H' term, since it contributes only to 
condensates, this leaves the terms multiplied by W, W and Ww The latter is in fact redundant, as we discuss below, 
so can be ignored. The key point is then that, in both the W and W terms, the lattice spacing appears in the 
combination 

Tr(i f E + E f i) . (18) 

If we now expand E about So in the following way, 

E = exp(iw T 3 /2) exp(i7? ■ f/f) exp(iw T 3 /2) , (19) 

then the quantity (|18fl . with A set to its final value d, is invariant under the spurionic symmetry 

n(x) -> -7r(x) , w -w , P -> -p, in <-> r- M =>• E E f , D M E ^ -D^E 1 , x' ^ • (20) 

It follows that terms even (odd) in pion fields must be even (odd) in loq. Since the only functions of luq which can 
appear are Cq and Sq, the former must be multiplied by an even function of pion fields, the latter by an odd. Thus at 
maximal twist (co = 0) the quantity (|18fl produces only odd powers of pion fields, and, in particular, has no vacuum 
expectation value. We also need to consider the combinations 

and Tr( X ' f E + E f x) , (21) 

which appear in the W and W terms. These are also invariant under the spurionic symmetry, but are independent of 
ujo, and thus must be even in pion fields. Now, when determining the consequences of the NLO terms, it is sufficient 
to expand about the LO vacuum, as we have done in 119|) . It then follows that the W and W terms only give rise, at 
maximal twist, to vertices involving odd numbers of pions. The physical vertices, which involve even numbers of pions, 
do not receive any corrections proportional to a and are thus automatically improved. Of course, the contribution 
linear in the pion field coming from the W term leads to an 0{a) tadpole, which can convert LO into NLO vertices 
with one less pion. However, since the LO vertices only have even numbers of pions, the resulting vertices with 0(a) 
corrections all contain an odd number of pions. Again, the physical vertices are not corrected until 0(a 2 ). 
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These arguments can also be extended to external sources, with the conclusion that physical matrix elements are 
automatically improved. Conversely, parity violating matrix elements are maximal at maximal twist. According to 
the general considerations of Ref. @, these results are expected for m q ~S> aAq CD . Our argument here shows that 
they hold also in the GSM regime in which m q ~ aAq CD . This does not require the use of an improved quark action, 
but it is essential to use a definition of maximal twist that sets ml ', rather than m, to zero. 

A possible concern is how these considerations generalize to extensions of the chiral Lagrangian incorporating other 
particles, e.g. baryons and heavy-light mesons. The point is that, although the Pauli and mass terms in the effective 
continuum Lagrangian © have the same chiral transformation properties, they are not proportional as operators. 
Thus, although they will always enter any extension of the chiral Lagrangian with the same form, their relative 
strength (Bq/Wq in the pionic sector discussed above) will depend on the quantity being considered. For example, 
the coefficient of a baryon mass term of the form Tr (BxB) and that of the corresponding term caused by discretization 
errors, Tr(BAB), will not be in the ratio Bq/Wq. Thus using the shifted mass ml will not, in general, remove the A 
term. It seems, then, that we should be concerned, when m ~ aAq CD , that the mass and discretization contributions 
are of the same size in general. The resolution of this concern is simply that the term linear in a is removed by 
working at maximal twist. In fact, we can see this in our calculations, because quantities such as J^, which do not 
vanish in the chiral limit, have a similar status in the chiral expansion as does the nucleon mass (although the details 
of the expansion differ). 

We conclude this section by addressing a few technical points. First, we discuss the W\q term. As noted above, this 
term is redundant: it can be transformed into a combination of the W, W and H 1 terms by the change of variables 

,5£ = ^ £ (Eit£-i). (22) 
J 

The shifted coefficients of the remaining three terms are W + W10/4, W + Wio/2 and H' — W\q. All physical quantities 
must depend on these combinations and not on Ww, W, W and H' separately. For this reason, we refer to them 
as "physical combinations" of LECs, despite the fact that they are introduced by discretization errors. Although we 
could remove W\q by this change of variables, we have kept it, both because we have found it to provide a useful 
diagnostic in the computations of matrix elements, and because it helps in understanding the impact of improving 
the underlying lattice theory, as we now discuss. 

We next consider how the effective Lagrangian changes if the underlying Wilson action, currents and densities 
are 0(a) improved, as this illuminates the nature of the LECs introduced by discretization errors. 4 The required 
improvement coefficients have been presented in Ref. 25|. As already noted, mass dependent 0(a) corrections are 
of higher order than we consider here. Thus the only improvement coefficients that are needed are csw, i-e. the 
coefficient of the lattice Pauli term in the quark action, and the standard improvement terms for the vector and axial 
currents proportional to cy and ca, respectively. We refer to Refs. [25|,|28| for details of these terms and discussion of 
how they can be determined non-perturbatively. If csw is set to its non-perturbative value, then the coefficient b\ in 
Eq. J2| vanishes, so that corrections to the quark effective action are of 0(a 2 ). It follows that Wo = 0, and the 0(a) 
shift in m c is absent, so that \' = X- The chiral effective Lagrangian then takes the same form as in Eq. (|I6[) . except 
that the factor of Wq in A should be dropped (with concomitant changes in the dimensions of various LECs), and, 
most importantly, W = W = 0. These latter two terms, which are proportional to a% and ap 2 respectively, are absent 
because there is no term linear in a in the underlying quark Lagrangian. The coefficient Wiq does not, however, 
vanish since it corresponds to 0{a) corrections to currents, and these have not yet been improved. Similarly, H' is 
non-zero. Finally, W is non vanishing since it is proportional to a 2 . If the axial current is also non-perturbatively 
improved, then one has in addition that Wio — (as can be seen from the result below for /a)- Improving the vector 
current by adding the cy term has no further impact on the LECs, since it turns out that cy does not contribute to 
the quantities we consider at NLO, as we discuss further below. In summary, if one non-perturbatively improves the 
action and currents, then the only remaining discretization errors are those of 0(a 2 ) proportional to W , as well as 
the H' term which contributes terms of 0(a) only to the condensate. Our subsequent results show that, in this case, 
while all physical quantities are then 0(a) improved, the unphysical parity-flavor violating matrix elements do have 
contributions at NLO which are proportional to a 2 /m q (see sec. IIII Hi below). 

Finally, we recall that the currents and densities can be obtained by taking appropriate functional derivatives of 



4 We thank the referee for suggesting that we consider the impact of improvement. 
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the action S — J d A xC x : 

*-K«? + B5)*- (23) 

S- = ~,S, S° = 4-„S, (25) 



is" 



where s = s° + J\=i sfc Tfc, an d similarly for p, and r M , although for the latter two we use only the flavor non-singlet 
parts. Note that these equations give the currents and densities in the twisted basis, which is that usually used in 
simulations. Our normalizations are such that the corresponding currents and densities at the quark level are obtained 
by taking the same functional derivatives as defined above of 

d 4 xC cS - ([r p + ty]/2 + 75 [r> - £ M ]/2) V + $(a + tj 5 p)^ (27) 

Note that the factors of 1/2 in eqs. <|23I26[1 imply that we are using Tk/2 in flavor non-singlet operators, while for the 
singlet we do not include the factor of one half. Thus, for example, P k — ^75^/2^ and S° = "ipip at the quark level. 

An important issue is the normalization of the currents. This can be determined, in the continuum and chiral 
limits, by enforcing appropriate Ward identities 29]. The currents we use at the quark level automatically satisfy 
these identities. This carries over to the currents in the chiral theory, since the matching maintains normalizations. 
The ratio of the scalar to pseudoscalar densities is also correctly normalized. The overall normalization of these 
densities, however, is a scheme dependent quantity, and this ambiguity is reflected in the chiral theory in the presence 
of the parameter Bg in \- Because of these considerations, the results we give below for matrix elements of these 
operators in the effective theory apply directly to lattice matrix elements, as long as the lattice operators include their 
matching factors. 

Away from the continuum limit, the explicit breaking of symmetries implies that one cannot, in general, normalize 
the currents and the ratio of densities in a universal way. Different choices of normalization condition will lead to 
results differing by 0(a), in general. We stress that this is true both on the lattice and in the effective chiral theory, 
as must be the case since the latter is supposed to represent the former. This is not, however, an impediment to 
our calculations. Since the effective chiral Lagrangian includes all terms that are allowed by the reduced symmetries 
of the lattice theory, our results contain, explicitly, all possible non-universal terms at the order we work. Thus in 
particular our results for matrix elements of currents and densities should hold for any choice of lattice operators 
which are correctly normalized in the continuum limit. 5 



III. CHIRAL PERTURBATION THEORY FOR GENERIC SMALL MASSES 



In this section we work out in detail the consequences of the effective chiral Lagrangian in the GSM regime. In 
fact, we are now in a position to state precisely how we define this regime. What we require is that 

> M' > a , (28) 

where M' is defined in Eq. (|13fl and A x = 4m f. This ensures that the LO terms in the chiral Lagrangian (|16fl dominate 
over the NLO terms. This condition can also be written using quark masses: 

Tfl 

1 » -7 — — > aA QCD . (29) 
Aqcd 



This argument fails, however, if one determines the normalization of currents and densities using the Schrodinger functional, as discussed 
in the cae of tmLQCD in Ref. 1251 . In that case there are additional O(a) contributions related to the boundary fields which are not 
accounted for by our analysis. It seems likely that our analysis could be generalized to include such boundary terms, but we have not 
worked this out. Thus we are assuming that the normalization constants are determined without the use of boundary fields, e.g. for 
Zy and Za using the method proposed in Ref. We thank the referee for pointing out this issue. 
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Note that our results for the GSM regime remain valid if a is smaller than M' . The converse does not, however, hold. 
In particular, if M' becomes as small as a 2 /A 2 then one enters the Aoki regime where NLO and LO contributions are 
comparable. 



A. The vacuum and Feynman rules 



The NLO terms cause a small realignment of the vacuum expectation value of E away from So- In fact, we know 
from continuum chiral perturbation theory that the Li terms in the potential do not realign the condensate. While 
this is possible in principle, given the physical values of the Li, it does not occur for quark masses in the range of 
interest to simulations. Thus the only realignment of the condensate is due to the W and W terms. 

We define the full NLO condensate to be 

(0|E|0)at LO = E m = exp(iw m r 3 ) = exp(i[w + e]r 3 ) , (30) 
so that e = u> m — loq. Minimizing the potential, we find 

16 



2 as (W + 2W'ac /M') , (31) 



The a 2 /M' term is not singular in the GSM regime, since M' is not allowed to become smaller than a. We note 
that, to the order we are working, we could make the replacements sq — > s m = sinw m and Co — * c m = cos uj m in the 
expression for e. 

Later on we will further replace So and Co, respectively, with the sine, s, and cosine, c, of a non-perturbatively 
defined twist angle to which differs from both ujq and Lo m by 0{a). This is allowed as long as we make the replacement 
only in NLO terms. Because of this, from now on we will use s and c without subscripts in NLO terms, and point 
out explicitly where the choice of twist angle is important. 

We now expand about the vacuum expectation value, defining physical pion fields by 

S = £ m Ep/j£ m , £ m = exp(iw m r 3 /2) , E p h = exp(«7? • f/f) . (32) 

We use an axial rotation, which has the symmetric form, since this is the rotation needed at the quark level in the 
continuum limit to undo the twist. It gives rise to pion fields that are in the physical basis, as we will see below. 
Inserting this expansion into the chiral Lagrangian (|16|) we find 

£ x = y(M' + AM')-f paVW 

+ ■ ^1 + ^M'L 45 + ^acW 

7T3 „ _ _ „ 16 „ ~ 7r 3 7T 2 eM' 
- y^TT • dpICjgasW + — — 

(M' + 4AM>) 16 , 

~ ^ ' 24/2 h n ^3 3^4 a s w + ■■■ K 66 ) 

where ir 2 = ff • if , and 

32 

AM' = j2 (M' 2 L 6S + acM'W + a 2 c 2 W') . (34) 

We stress again that s and c could equally well be s$ and cq, or s m and c m , at NLO accuracy. With one exception, 
we show explicitly in 133|l only terms which we need for our calculations below. In particular, we have not included 
vertices proportional to L\ and L2, nor four pion vertices with two derivatives. These are unchanged from continuum 
chiral perturbation theory [the rotation in Eq. (|32l) cancelling in these terms] , and so lead to NLO contributions of the 
same form as in the continuum. The exception are the vertices involving four pions, which we do not need explicitly 
below, but which we included to show the type of flavor breaking which occurs. 

These results illustrate the generic features discussed in the previous section. First, parity and flavor conserving 
terms of 0(a) are also proportional to c and thus vanish at maximal twist. Second, the 0(a) terms which do not 
vanish at maximal twist (i.e. the three pion vertices proportional to as) violate parity and flavor. The factor of s 
ensures that they vanish when /1 — > (when the flavor symmetry is restored). Finally, flavor breaking, but parity 
conserving contributions (i.e. the additional mass-term for 7r 3 and the 7r 2 7r 2 vertex) are proportional to a 2 and to 
s 2 oc {J 2 . All these results are expected on general grounds for M' ^> a, but also hold, as we see here, for M' ~ a. 
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B. Currents and densities in the twisted basis 



At LO, the currents and densities take the usual form (although other authors use different normalizations): 

K,lo = yTr (^[E^E + EfyEt]) , (35 ) 

<lo = Y Tr ~ E9 ^ f ]) , (36) 

S£ o = -^Tr(T fe [E + £t])=0, (37) 

S£ = -^Tr(£ + £t) , (38) 

P|o = ^Tr(r fc [E-St]) , (39) 

= ^Tr([E-£t]) =0 , (40) 

where fc = 1, 2, 3. The vanishing of the isovector scalar and isoscalar pseudoscalar densities is a property of SU(2), 
and holds also at NLO, so we do not consider these quantities further. 
At NLO, the vector and axial currents become 



V£ = V* LO (l +C) + L U L 2 , L 9 terms , (41) 
4aWn 



K = < i0 (l + C) + ^^a„P* + L X ,L 2 , L 9 terms , (42) 



where k = 1, 2, 3 and 



C = ~^2~ Tr ^ S + S X ] - 2Bq j 4 Slo ■ (43) 

We do not give the form of the L\, L 2 and Lg terms since these are unchanged from the continuum p7T |. 

For the scalar and pseudoscalar densities at NLO, if we drop contributions proportional to the sources s or p, which 
give rise only to contact terms in correlation functions, we find 

S° = S" LO (l +V) - 8B (cM'H 2 + aH') , (44) 

P a = P£ (l + V), (45) 

P 3 = P'loi 1 + + 4iB sM'H 2 , (46) 



where a = 1, 2, and 

V = -^Tr[^ED,Et] + ^Tr( X 't S + ^ _ i|^ . (47) 



Note that the terms in the chiral Lagrangian proportional to L 5 , Lj and L% give no contribution to these densities, 
due to the properties of SU{2) matrices. The constant terms proportional to M' and a in 5° and P 3 contribute only 
their vacuum expectation values. 

We stress that all expressions in this subsection so far are written in terms of E, whereas to use them we need to 
change variables to Ep^. To do this we first need the results: 

ya _ „ I/O _ ,3afc 2b t/3 _ t/3 

v \x,LO — °m V^LO e s m A p,,LO > v n,LO ~ v n,LO J 

A a — r A a — e 3ah v V h A 3 — A 3 

Sip = C mSlo ~ iSm^pLO ' PlO ~ c mP\o ~ * s m'S'" /2 , P^q = P£q , (48) 

where a, b = 1, 2, and the "hatted" currents and densities on the right hand sides of these relations take the same 
form as the LO currents and densities in eqs. Ij35l4t)|l but with E replaced by Epfc. For example 

Klo = Y Tr ( r2 PphWU + Kh 9 ^P»}) ■ ( 49 ) 
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These are LO currents and densities of the chiral effective theory expressed in terms of physical pion fields. The 
transformations in l|48ll are identical to those at the quark level between the physical and twisted basis. We stress 
that to obtain these results it is essential to relate E to T, p h using the axial transformation of Eq. 1(32(1 . and not, say, 
a left- or right-handed transformation. 
Finally, we will also need the results 

Tr[ X 't S + EV] = M'Tr[Ep/j + Ejj + 0(M'e) = + O(M'e) (50) 

Tr[D M E£»^Et] = Tr^E^Ejj , (51) 

which allow us to express C and T> in terms of physical pion fields (with the spurion \' set to its final value). We can 
drop the O(M'e) term in the first line as this contributes only at NNLO. 



C. Defining the twist angle 

Using the expressions for the currents and densities, we can determine the twist angle non-perturbatively. In the 
continuum limit, if the input twist angle is to — tan (///m), the physical currents and densities are given by 

V^cV^ + e^sAl, V 3 = V 3 , 

S a = cS° + is2P 3 , P 3 = cP 3 + isS°/2 , P a = P a , (52) 



where a, b = 1, 2, c = coslu and s — sinuj. As already noted, these have the form of the inverse of the transformations 
of the LO operators, Eq. (|48|l . except that here the twist angle is to rather than to m . Note that we can unambiguously 
create charged pions using P a and neutral pions using A , since these operators are invariant. 

In the continuum, we can determine to by the condition that there is no flavor or parity breaking in the physical 
basis. Different ways of enforcing this all lead to the same result for to. On the lattice, however, discretization errors 
imply that the different definitions will lead to results differing by 0(a). 

We will take as our canonical definition of to that obtained by enforcing 

(V 2 (x)P 1 (y))=0. (53) 

The result for to depends, at 0(a), on the distance \x — y\. We enforce the relation at long distances, where the 
single-pion contribution dominates. Using the definitions in l(52[). we can express the condition as 

(v?(*)PHy)) ( , A , 

tanw = , , . . . — , . ,. . (54) 
(Ai(x)Pi(y)) V ' 

We stress that, at long distances, this gives the same result as that involving a divergence, 

(d^(x)P\y)) 

tanw = — — /! ; r — t——- . (55) 

(d^(x)pi(y)) v ' 

since the factors of <9 M acting on the pion propagator cancel. Both definitions l|54|l and l|55|l have been used in 
simulations (see, e.g., Refs. [SME1)- 

To evaluate (|54() at LO we use the results in Eq. 148(1 . At this order, only A 1 ^ LO and P 1 — P 1 couple to the 
single pion state, and we find that lu = Lu m . Since tu m — luq at LO, we also have that lu = luq. This shows that 
the non-perturbative definition l(54|l automatically includes the shift from m to ml discussed above. In particular, if 
lu = 7r/2, then m! = up to corrections of 0(a 2 ). 

The NLO calculation of lu is simplified by the fact that most contributions cancel in the ratio (|54|l and do not 
change the LO result. This holds for the factors of 1 + C in eqs. ((41142(1 (evaluated with = 1), as well as one- loop 
wave- function renormalization, and the one-loop corrections to the coupling of A 2 ^ LO and P 1 to the pion. The L\, 
Li and Lg terms do not contribute at this order. The sole non-trivial contribution is that from the TUio term in 
and we find 

s m 4aW 10 
tanw= — , d = — -g— , (56) 

Cm+O J 
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This result can be rewritten at NLO in a number of useful ways: 

s = s m - 5s m c m , c = c m + Ss^ , lo - u m — ~5s m . (57) 

Note that in each of these relations one can substitute s for s m , etc., in the sub- leading term on the right hand side. 

As discussed earlier, Wio cannot appear alone in a physical quantity. This is not a concern here, however, since 
uj m is not physical, i.e. not directly accessible through a non-perturbative calculation. What must be physical is 
the difference between results from alternative non-perturbative definitions of lo. In particular, we expect loq to be 
physical since it contains the physical quark masses. Indeed, the difference 

uj - oj = (W + Wio/4 + 2acW'/M') (58) 

is given by a physical combination of LECs. 

Another non-perturbative definition of lo can be obtained by enforcing the continuum relation 



(V?{x)Al(y)) = 0. (59) 



This leads to the result 



taW= (WFpp- (60) 

In lattice simulations this relation has been used along with Eq. (|54|l in order to determine u> without knowing the 
normalization of the vector and axial currents Here we do know these normalizations, so this relation gives, in 
principle, an independent determination. We find, however, that the result for lo is identical to that obtained from 
B54JI . In fact, this could have been seen in advance since, when using the single pion contributions, enforcing both 
|[5DJ and (3D amounts to requiring that (OIV^tt 1 ) = 0. 

We do obtain a different non-perturbative result for the twist angle at non-zero lattice spacing if we enforce 
(S°(x)Al(y)) = 0. This leads to 

i(S°(x)Al(y)) 

tanWpE ^W' (61) 
where we have denoted the new angle lop (and will use cp = coswp, etc.). A straightforward calculation leads to 

4as(AW + W w ) 

uj p -uj= — 75 (62) 

at NLO. The difference between u> and ujp contains a physical combination of LECs, as it must since they both can 
be calculated non-perturbatively. 

Another way of stating this result is to use uj obtained from l|54|l and calculate the coupling of the physical scalar 
density to 773. The mismatch of twist angles implies a non-zero result, which can be expressed as 

(S°(x)Al{y)) = ~8ias(4W + W w ) 
(P 3 (x)Al(y)) P 

where again we take only the single-pion contribution. The result would, of course, have vanished if the physical scalar 
density had been constructed using ujp rather than lo. We note that (|B3|l [or (|BT|l ] provides a method of calculating 
one combination of the LECs that are introduced by discretization errors. We return to this point in more detail 
below. 

In summary, using the currents and densities we have obtained two non-perturbative definitions of the twist angle. 
These are equivalent in the continuum limit, but differ by 0(a) away from this limit. This 0(a) ambiguity in the twist 
angle is maximal at maximal twist. Nevertheless, as we show below, it affects masses and physical matrix elements 
only at 0(a 2 ), and thus does not contradict automatic 0(a) improvement. Only flavor-parity violating quantities are 
affected at 0(a). 
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D. Constructing the physical currents and densities 



Using the NLO result for uj denned by l|54|l . which we take as canonical from now on, we can now explicitly construct 
the physical currents and densities, using the relations in Eq. 152|) . We find for the currents (with a, b = 1, 2): 



b 

/j,,LO 



AasW, 



10 Jiab 



V: = V; 



(1 



KM 1 



KM 1 



c)- 
-c) 



p 

L\, L/2, Lg terms , 
AasW w * ab 



P 

2B P a " 



B P 

AacWy 
^P 



e MO daPlo+Li,L 2 ,L 9 terms, 



Wo 



L\, L2, Lg terms , 



L\, L27 Lg terms , 



where C is the same as above, but we now write it using physical fields: 

16 



2B p 



-(M'L 45 + acW)S° LO + asW2iPl 



(64) 
(65) 
(66) 

(67) 
(68) 



The L\, L2, and Lg terms take exactly the same form as in continuum chiral perturbation theory, but now expressed 
in terms of This is because they rotate exactly like the LO currents under the axial rotation. Since we are 

interested in flavor-parity breaking contributions we do not give these terms explicitly. 
For the densities, the results are 



S° 



S° 



LO 



(1 



P a = P a LO {\ 



P 3 = P| (1 



■V)- 

■T>), 

■V)- 



AasWi 
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P 



4iB asH' , 



4asW 10 iS° LO 



where we now express T> in terms of physical fields: 

4£ 45r 



V 



P 



-Ti[D^ ph D^l 



P 
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-SB (M'H 2 + acH'), 



phi 



2B p 



-(2M'L 6S + acW)Sl + asW2iPl a 



(69) 
(70) 
(71) 

(72) 



For both currents and densities, the mismatch between u> and uj m leads to contributions on the right hand sides 
proportional to W10 . Note that the currents and densities themselves need not be composed of physical combinations 
of LECs, since they are not directly measurable. It is their matrix elements, such as those we compute below, which 
are physical. 

One test of the results given above is that the single pion matrix elements of should vanish, i.e. the physical 

vector currents do not couple to single pions. There are contributions both from the A h ^ LO and d^P^o terms and 
these do cancel. In effect, our condition for determining u> has enforced this result. 

There is a subtlety concerning the choice of twist angle. The results above for the physical currents and densities 
hold as written for our canonical choice, to. If instead we had used ujp to define the rotation between twisted and 
physical bases, then some of the explicitly parity-flavor breaking terms are changed. In particular, in the terms 
proportional to A* in Eq. to V* in Eq. l(So]l. to P 3 in Eq. l(o§)l. and to S° in Eq. (JJTJ, W 10 would be replaced 
by —AW . The other terms proportional to W\q are not, however, changed. In our subsequent results, these changes 
would only impact the flavor-parity violating matrix elements. 

Finally, we note the same general features in these results as observed earlier in the Feynman rules. Discretization 
errors which do not violate flavor and parity come with factors of c, and thus vanish at maximal twist. Flavor-parity 
breaking terms, however, are proportional to as. 



E. Pion masses 



With the Feynman rules, currents and densities in hand, we now turn to the predictions for masses and matrix 
elements. We begin with the charged pion mass, which we find at NLO to be given by 

ml 12 = M' + i| (m ,2 (2L 68 - L 45 ) + M'ac(2W -W) + 2a 2 c 2 W^ + 1-loop . (73) 
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Here the one-loop contribution is unchanged from that in the continuum as long as the result is expressed in terms of 
the LO pion mass-squared M' . This is because it involves only LO vertices and masses, which themselves have the 
same form as in the continuum. To be completely clear on this point we quote the result in this case 

M'\ ( M' 



where is the renormalization scale. Since our emphasis is on discretization errors, however, we will not give 
explicit expressions for the 1-loop contributions to other quantities below. They can be found in the original work on 
two-flavor ChPT [27j . as well as in more recent works including the extension to the partially quenched theory, e.g. 

Refs. muni. 

There are various checks on our result l|73l) . When s — ► 0, it agrees with that of Ref. 0], where m 2 was calculated 
to NLO in the untwisted Wilson theory. It goes over to that of Ref. 0] if we take a ~ M' 2 . And, after considerable 
algebra, it can be shown to agree with the result of Ref. [l4| , where the same quantity was obtained. 

Various features of the result are noteworthy. First, it depends only on physical combinations of LECs, as required. 
Second, it is automatically 0{a) improved at maximal twist (c = 0), or after mass averaging. Finally, the a 2 correction 
provides a further shift to the critical mass. This shift vanishes, however, at maximal twist. Thus we find the somewhat 
surprising result that, at maximal twist, the charged pion mass differs from the continuum only by terms of NNLO 
in our expansion. This is not to say that the result is 0(a 2 ) improved, because there are contributions proportional 
to M'a 2 . 

The flavor breaking in the pion masses at NLO is given solely by the analytic contribution, and we find (in agreement 
with Ref. H3) that 

,2 „,2 _ ^2 2 2 



7 ' " ~ f 2 



= --pra 2 s 2 W (75) 



J* (m + ay + jx* 

To obtain the final form we have used the fact that we can replace s with s a t the order we are working. We know 
on general grounds that this splitting must vanish quadratically in a/x, because it does not violate parity. Naively, 
then, one might have expected the splitting to arise first at fourth order in our expansion (NNNLO). What our result 
shows is that, in fact, there is mass dependence in the numerator such that the effect is of NLO. 

Calculating the pion mass splitting in practice is complicated by the fact that the neutral pion propagator includes 
disconnected quark contractions. It would nevertheless be an interesting quantity to determine. For one thing, it 
would give an indication of the size of discretization errors. Furthermore, as noted in Ref. |14| . determining W 
gives information about the phase diagram in Aoki region. In particular, if W < there is an Aoki phase on the 
untwisted Wilson axis , w hile if W > there is a first order transition on the Wilson axis extending out into the 
twisted plane [Hjl IbH Il5l Il9j . It is perhaps surprising that a simulation in the GSM regime can give information 
about the Aoki regime. The reason for this result is that the same terms in the chiral Lagrangian are responsible 
both for the mass splitting and for determining the phase structure. In fact, the first form of the result, 1)75(1 . holds 
also in the Aoki regime, although the final form, (|76|) . does not. We discuss this further in the next section. 

F. Decay constants and PCAC quark mass 

Using the expressions for the axial current and pseudoscalar density in the physical basis, we can determine their 
one-pion matrix elements up to NLO: 

Sa = f jl + -^[2M'L 45 + ac{2W + W w )] + 1-loop J , (77) 

fp = fB jl + js [M'(4L 68 - L 45 ) + ac(2W - W)] + 1-loop j . (78) 

Here, as for the pion masses, the one-loop term is the same as in continuum chiral perturbation theory |2?| and we do 
not give it explicitly. At this order there is no flavor breaking — this enters first at 0(a 2 ), which is NNLO for the decay 
constants. The result shows the expected automatic 0{a) improvement at maximal twist or under mass averaging, 
and the appearance of a physical combination of LECs. 
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In fact, the quantities that we have discussed so far — w, m%, Ja and fp — are not independent. To see this we first 
consider the so-called "PCAC quark mass" defined through 

_ (d^(x)P a (y)) 
mpCAC = 2(P°{x)P*{y)) ' (79) 

where a = 1,2. Note that it is not the physical axial current which appears, but rather that in the twisted basis. 
This quantity is of interest because in the continuum limit it gives the untwisted component of the quark mass. We 
evaluate the correlators for \x — y\ — > oo so that the single pion contribution dominates. At LO we find 

™ L p°cac = m'= , (80) 

tan ujq 

showing that this quantity automatically includes the 0(a) offset in the untwisted mass.. At NLO we find 

„ fNLO( m 2 \NLO 

C JA \ m TT a ) / Q1 \ 

m PCAC = yWTo > ( 81 ) 

where the superscripts indicate that the full NLO expressions must be used. In this result, we cannot change the 
overall factor of c = cosw into, say, c m or cq. 

The result (|81|l also shows that mpcAC is not independent of the quantities considered so far. In particular, using 
the condition mpcAC = to determine maximal twist is equivalent to setting tanw = oo in 154(1 . Both rely on the 
vanishing of the coupling of the axial current in the twisted basis to the pion. 

To proceed further we recall that the lattice symmetry has an exact "PCVC" relation 

d;V^(x) = -2fi e 3kl P l (x), M = l,2,3,. (82) 
where 5* is the backward lattice derivative, P k the local bare pseudoscalar density, and V the point-split current: 

■y k k 

V£(x) = ^(z)^ - 1) T —U{x 7 ^{x + a/2) + $(x + aA)(7* + 1) \u{x, • (83) 

The same relation must hold in the effective quark and meson theories, to all orders in our expansion, if we normalize 
our operators to maintain continuum Ward identities. In particular, at the meson level, we have 

d^(x) = -2ne 3kl P l (x). (84) 

Using this, and the definition of the twist angle l)54[l. we find 

_ {d,Al{x)PHv)) (d^(x)PHy)) n 
mpcAC {d ^ V 2 [x)P i {y)) x 2 {P\x)P^y)) tanw ' 

This relation should hold to all orders in tm^PT, and we have checked that it is valid at NLO. Using Eq. 1(81(1 . we 
see that /a, fp, m% and u are related. 



G. Parity conserving matrix elements 



We next consider other flavor-parity conserving matrix elements that can be calculated in numerical simulations. 
The first example is pion vector form factor, which can be obtained from the matrix elements 

Mpi)|U> m ( P2 )) . (86) 

We find that, at NLO, the result is unchanged from that in the continuum, which is given, for example, in Ref. |27j . 
This holds true at any twist angle. It follows because the contribution from the factor C in 1(68(1 cancels that from 
wave-function renormalization. 

It is perhaps surprising that there are no 0{a) terms even for untwisted Wilson fermions. We know in this case 
that to remove 0(a) terms in general we need to add to the current a term containing the tensor bilinear. One can 
show, however, that this contributes to the pion matrix element terms suppressed by 0(aq 2 ) compared to the LO 
contribution. While of O(a), these are NNLO corrections in our power counting scheme. 
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Our next quantity is the scalar form factor of the pions. 
from the matrix elements: 



This are flavor conserving at NLO, and can be obtained 



(7T k { Pl )\S°\TT k (p 2 )) = 2B 1 1 + [ - q 2 L 45 + 4M'{2L 6S - L 45 ) + 2ac(2W - W)\ + 1-loop 



(87) 



As usual, the one- loop contribution is the same as in the continuum |27| . Unlike the vector form factor, here 0(a) 
improvement occurs only at maximal twist. 

We consider here also vacuum expectation values of 2iP 3 and S°. The former is the most interesting since it is 
calculable on the lattice and, at maximal twist in the continuum limit, it gives the physical condensate (0\uu + dd\0). 
It is calculable because P 3 does not mix with the identity operator in the chiral limit. Away from this limit, the 
mixing is proportional to /i/a 2 and can, in principle, be extrapolated away. By contrast, 5° mixes with the identity 
operator for all quark masses, and thus is very hard to calculate on the lattice. Nevertheless, we quote results for it 
as they illustrate an interesting theoretical point. 

Our results for these condensates at NLO are: 



(2iP 3 ) = -2fB s \l + — [M'(8L 68 + Hi) + ac(AW + Ww)] + 1-loop 



(S°) = -2f 2 B Q c \ 1 + [M'(8L 68 + H 2 ) + ac(4W + W w )] + 1-loop \ + 8B Q a(W 1Q - H') 



(88) 



(89) 



Here the relative 1-loop correction is of the same form for both quantities and of the same as for the continuum 
condensate. We stress that the inclusion of the H' term in the chiral Lagrangian is essential to make the result for 
(S°) physical. We note that (P 3 ) is O(a) improved at maximal twist, which is expected because it is proportional 
to the physical condensate (in the chiral limit). We note also that the M' term in (P 3 ), since it is multiplied by the 
overall factor of s, is indeed proportional to fi. Thus H 2 is the constant which contains the quadratic divergence. On 
the other hand, S° is not a multiplicatively renormalizable operator, as it mixes with the identity operator, even in 
the chiral limit. Thus the general considerations of Ref. [5j do not require that it be automatically 0(a) improved, 
and indeed it is not, as shown by the last term. 



H. Parity violating matrix elements 



We finally consider simple examples of unphysical, parity violating quantities which vanish in the continuum due 
to flavor and parity symmetries, but are present in tmLQCD at O(a) since they are not automatically improved. The 
first such quantities are axial form factors of the pions. These are obtained from the matrix elements 



(Taluks), kal^Ka), folios), 



(90) 



where a = 1,2. Calculating these in lattice simulations is straightforward in principle as they involve only single 
particle states. The main complication is that there are quark-disconnected contractions in addition to the usual 
quark-connected contractions. 

At NLO, the contributions to these matrix elements are of two types, as shown in Fig. (JTJ: the direct two pion 
terms in the axial current (both the A^qP^q and the V^o terms), and a pole term in which the axial current creates 
a single pion which then connects with the two external pions through the three pion vertices in C x . There is no 
contribution from the L\, L 2 , or Lg terms since these do not break flavor or parity at this order. We also do not 
need to include wavefunction renormalization since the leading contribution is of NLO. There are also no one-loop 
contributions. Our results are (with a = 1, 2): 
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(a) Direct term 



(b) Pole term 



FIG. 1: Generic types of Feynman diagrams that contribute to the two pion matrix elements of the physical axial current. 
The arrows label only the momentum flow. The direct term comes from the two pion vertices in the axial current. The pole 
term arises from contracting the single pion term in the axial current with the three pion vertices in C x . 

where q = pi — Pi- A check on the results is that all quantities are composed of physical combinations of LECs. Note 
that the overall factor is maximized at maximal twist. 

The final quantities we consider are the pseudoscalar form factors of the pions. The required matrix elements are 
similar to those for the axial currents, and have both direct and pole contributions. Our results are (with a = 1, 2): 
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(94) 
(95) 
(96) 



These are very closely related to the corresponding axial matrix elements. 

The results for the parity-flavor violating form factors do depend on the choice of definition of twist angle. The 
results above are for our canonical choice, us. If instead we use usp to rotate to the physical basis, then we find 
that in Eqs. JHH, EEl and (JHEJ [though not in Eqs. and ffity] the factors of W w /A are replaced by -W. In 

particular, this means that all the matrix elements of the physical pseudoscalar density do not depend on W\o with 
this definition of twist angle. This explains a puzzle that appears when considering the application of our results 
to an underlying lattice theory in which the quark action is non-perturbatively improved, but the currents are not. 
Since the pseudoscalar density does not require improvement through NLO in our expansion, one would expect that 
its parity violating matrix elements would have no terms linear in a. As discussed above, however, improving the 
quark action implies setting W = W = but leaving Wio non-zero, and so the results of Eqs. (|95|l and l)96|) do not 
vanish. The puzzle is explained by noting that the twist angle is determined using the axial and vector currents, 
which are unimproved. If instead one uses the pseudoscalar and scalar densities, i.e. uses asp, then, as just noted, the 
W w dependence is absent as expected. 

The matrix elements of this subsection provide a way of determining the LECs associated with the discretization. 
Using the momentum dependence one can, in principle, separately determine the physical combinations W — 2W, W' 
and W + Wio/4. Furthermore, since only two of the expressions in square brackets in (|91I93|) are independent, and 
similarly for (|94I96(I . there are a number of predictions implicit in these results. We stress that these predictions hold 
only at NLO in tm^PT, and will be broken at higher orders. 



I. Summary of results in GSM regime 



In this final subsection we summarize what we have learnt about the GSM regime using tm^PT, and describe one 
way in which these results could be used in simulations. We do not discuss the technical problems that can arise 
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when calculating the quantities we have discussed using lattice methods (e.g. the difficulty in simulating disconnected 
contractions). We simply assume that these difficulties can be overcome. 

We imagine that a simulation is done with particular values of fj,o and too, and that all the needed Z factors have 
been determined. Thus, given Zp, we know the renormalized twisted mass, fi. We do not, however, yet know the 
shifted untwisted mass to'; for this we need to determine the (shifted) critical mass with 0(a) accuracy. Our only 
assumption about to' is that it (and /i) are such that we are in the GSM regime. Our aim is to determine as many as 
possible of the parameters in tm^PT, to make predictions about physical quantities, and to provide consistency tests 
of tmxPT. 

One possible way of proceeding is as follows. 

1. We determine u> non-perturbatively using Eq. I|54|l. This determines the untwisted quark mass including the 
0(a) shift using 

' '" - -[I + 0(a)}. (97) 



tan wo tana; 

The relative uncertainty of 0(a) (which corresponds to an absolute uncertainty of 0(a 2 ) in the GSM regime) 
means, however, that the Aoki region cannot be resolved — it collapses to a point at this accuracy. 

2. To obtain to' to greater accuracy we next calculate top using Eq. 1|) . From the difference top—uj, or equivalently 
from the ratio given in Eq. (|63[) . we can determine a(4VF + Wio)/ f 2 . This is interesting in its own right, as it 
not only gives an indication of the size of discretization errors, but it also allows us to partially "correct" for 
the difference between cu and ujo- In fact, for maximal twist (uj = ±7r/2) we have, using Eqs. and 

ujq = oj P [l + 0(a 2 )] (Maximal twist) . (98) 

In this case we now know to' with a relative uncertainty of 0(a 2 ), and an absolute uncertainty of 0(a 3 ). We 
stress that this accuracy is to be understood within our power counting scheme, in which, for example, shifts in 
to' proportional to a/x 2 , which are expected to be present, are of 0(a 3 ) and so are too small to be included. 

Because of this result, one might argue that it makes more sense to use wp as the canonical choice rather than 
lu. One reason not to do so is that the calculation of uip is more difficult as it involves disconnected quark 
contractions. 

3. To obtain u>q (and thus to') at NLO accuracy for arbitrary twist angle we next calculate the pion mass splitting. 
This indicates the size of a 2 corrections, and specifically determines the combination a 2 W / f 2 . We can then 
obtain ujq using 

m l 3 ~ m l a rnn , 
ujq=ujp §- a. _ (99) 

toj tan u> 

Here a = 1, 2, and the mass of any of the pions can be used for m n in the denominator. 

At this point, we know where our simulation lies in the to', fi plane with an 0(a 2 ) relative uncertainty. Thus we 
know the physical quark mass m q = [i/sq = M' /(2Bq) with similar accuracy. This resolution is good enough to 
resolve the Aoki region (ml ~ /j, ~ a 2 Aq CD ), so we know what parameters to use if we want to study that region. 
Indeed, as discussed above, knowing the sign of W we can predict the phase structure in the Aoki region. 

4. We now use one of the parity-violating matrix elements to determine the final, linearly independent, unknown 
combination of LECs associated with discretization errors, namely 2W — W. For example, we could use 

2i(ir a (p 1 )\P 3 \Tr a (j }l )) 16as(2W-W) , . . 

(7r a 0i)|5°|7r°(pi)) r 

Here we have taken p2 = pi so that q 2 = 0, and normalized with the scalar matrix element to remove the overall 
factor of B . 

Given this result, we can make absolutely normalized predictions for the ratios of all the other parity violating 
matrix elements of P k to the pion matrix element of S° , including their q 2 dependence. Similarly we can predict, 
including normalization, the parity violating axial form factors given in Eqs. (|91I93|) . These predictions provide 
a test of the applicability of tm^PT at NLO, as well as an indication of the size of discretization errors. 
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5. Finally, we calculate the physical quantities m 2 , /a, fp, (P3) and (ir\S,V\ir). Here we distinguish between 
maximal and non-maximal twist. 

(a) Maximal twist. All these quantities are automatically O(a) improved, i.e. they have exactly the same 
dependence on m q as in the continuum up to 0(a 2 ) corrections, within our power counting scheme. Our 
results show that this automatic improvement holds even if one is working in a region where m q ~ aAq CD . 

We stress that this automatic 0(a) improvement holds irrespective of whether we use to or up or some 
other choice for the twist angle differing from lu by 0(a). This is because an O(a) uncertainty in lu leads 
to an uncertainty in physical quantities of size as, so that they are still 0(a) improved. 

(b) Non-maximal twist. Here we find that, having determined the LECs, we can completely remove 0(a) 
errors in the physical quantities that we have calculated. These errors can simply be subtracted or divided 
out. For example, the quantity 

l + l6ac{2W ~W)/f 2 

is predicted to have the same dependence on m q as does m 2 , in continuum %PT, up to 0(a 2 ) errors. 
This is an amusing, though perhaps academic result. In practice, when using tmLQCD, one should clearly 
work at maximal twist and avoid the need for such corrections. Nevertheless, it might be worthwhile using 
tmLQCD to determine the LECs (which, for small enough quark masses, are mass independent, though 
they do depend on the lattice spacing and gauge action) and then try to correct results previously obtained 
with unimproved untwisted Wilson fermions. 



IV. THE AOKI REGIME 



As noted in the introduction, it is of practical interest to study how the results obtained in the GSM regime 
change when one enters the Aoki regime. We recall that the latter is defined at the quark level as the region where 
m q ~ a 2 A QCD , or in tmxPT by M' ~ a 2 /A 2 . 

In the Aoki regime, competition between terms of size M' and a 2 / A 2 , leads to a non-trivial phase structure, with 
first order transition lines ending at second-order endpoints. This has been discussed extensively in Refs. [Tj, [Tj, Il3 
and we do not recapitulate the analysis here. We will present results which hold throughout the Aoki regime except 
on the phase transition lines themselves. 

Before presenting results we would like to make a general comment on the concept of 0(a) improvement in the 
Aoki regime. The improvement program is predicated on the assumption of small discretization errors. In particular 
the discretization effects are assumed not to cause large changes in the vacuum. This assumption clearly breaks down 
in the Aoki regime, where, in general, the direction of the quark condensate is determined by a competition between 
quark masses and discretization errors and the result differs by angles of 0(1) from the continuum theory. Thus in 
taking the limit a — > 0, quantities do not change in a manner which is perturbative in a, and the language of 0(a) 
improvement cannot be applied. This point has been emphasized in Ref. |3|]. 

A possible exception to this discussion has been raised in Ref. [21j . They argue that O(a) improvement at maximal 
twist remains valid in the Aoki regime, as long as an appropriate choice of twist angle is used. We return to this point 
at the end of this section. 

Irrespective of this general question, one can simulate in the Aoki regime and it is useful to obtain the predictions of 
tm%PT. We stress that, even though the language of 0(a) improvement may not apply, one can still use Symanzik's 
continuum effective Lagrangian to study the theory in the Aoki regime. 

Now we turn to the results. We will keep the discussion brief, minimizing technical details, since much of the work 
is a straightforward generalization of that in the GSM regime. 

As explained in Ref. |15|. we do not need to augment the chiral Lagrangian used in the GSM regime in order to 
study the Aoki regime. Instead, since the power counting is now p 2 ~ m q ~ a 2 Aq CD , and we work only at LO in this 
expansion, we can drop some of the NLO terms used above. In particular, to study the vacuum alignment and pion 
masses we need only keep the W term — those proportional to W, W and the Li can be dropped. When we consider 
the currents and densities, however, a new feature emerges. Since the functional derivatives defining these quantities 
effectively "eat up" one power of p 2 or M', terms which were of too high order for pion masses (such as the W term 
of size aM' ~ a 3 ) now need to be kept. Indeed, for each quantity that we consider, a careful study is necessary to 
determine which terms coming from the GSM regime can be consistently kept. 
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As already noted, the fundamental difference between the GSM and Aoki regimes is that the condensate is no longer 
closely aligned with the quark mass. In other words, e = uj m — uj is no longer of 0(a), but instead is generically of 
O(l). An indication of this is that the W contribution to e in Ea. 1)31(1 is proportional to a 2 / (f 2 M') and thus of 0(1) 
in the Aoki regime. The precise alignment of the condensate is determined by a quartic equation which has been 
given in various forms in Refs. |13|, [lj, lla. l21| . It can be written: 

M' S m(u m - lu ) = M' sine = - 32 ^ 3 ^ W ' + 3 

f 

Note that one solution is uj m — ujq — tt/2, so that e = at maximal twist. The corrections to this result are of 0(a), 
but we do not control these since they include contributions from a 3 terms in the chiral Lagrangian. 

We now work our way through the quantities considered in the previous section, concentrating on how they differ in 
the Aoki regime. The expressions for the currents and densities, 1(4214711 . are unchanged, except that the contributions 
proportional to M' or p 2 can be dropped. We do, however, control the 0(a) corrections proportional to W, W and 
Wio- The change of variables in Eqs. 14811 and ((50(1 also goes through unchanged, but that in l(51|) is replaced by 

Tr[x n S + £V] = Af'coseTrpV, + E+J + iM' sin eTr[n(E ph - sjj] . (103) 

In fact, this change only matters if this term contributes at LO, and thus only affects the pion masses, which are 
discussed below. For all other quantities the M' contributions are of higher order than we control in the Aoki regime. 

Because of these considerations, the determination of the twist angles to and up goes through unchanged. They 
differ from each other, and from uj m , by 0(a), with the differences given in Eqs. 1(57(1 and ((62() . The expressions for 
the physical currents and densities, Eqs. 1(64172(1 . thus remain valid, except that the terms proportional to M' and 
p 2 can be dropped. It then follows that the GSM expressions for /a, fp, the parity conserving matrix elements and 
the condensates remain valid in the Aoki regime, except that we should drop the M' and 1-loop contributions. We 
collect the main results here for clarity 

/A = / jl + -^ac(2W + Wio) + 0(a 2 )| , (104) 
fp = fBo jl + ^ac{2W -W) + 0(a 2 )| , (105) 
(n k ( Pl )\S°\7r k (p 2 )) = 2B Q jl + -^2ac{2W — W) + 0(a 2 ) J , (106) 



(0|2iP J |0) = -2f z B s jl + — ac(AW + W w ) + 0(a 2 ) j . (107) 

As in the GSM regime, we can replace s and c with sp and cp, respectively, to the accuracy we work. The vector form 
factor, which we do not display, is simply given by the leading order term in continuum xPT with no 0(a) correction. 
The results for the pion masses are changed more substantially. These have already been calculated in Refs. 0, 
so we only comment on the relation to our results in the GSM regime. In the Aoki regime, we need only keep 
the LO mass term and the W term from 1(16(1 . For the former, Eq. 1103(1 shows that the contribution to m 2 changes 
from M' to M' cos e in the Aoki regime. Combining this with the W contribution, which is unchanged from the GSM 
regime, we find (in agreement with Refs. |l4 Il5l l2llp 



in 



32a 2 c 2 W 3 , p, 



M' cos e + — — tj— — h 0(a 3 ) = — + 0(a 3 ) . (108) 



We have used Eq. I(102|) to obtain the second form. 

The mass-squared splitting comes only from W term in the chiral Lagrangian. It is easy to see that its contribution 
takes the same form as in Eq. ((3311 . with s being s m . But since s and s m differ only by 0(a), one can use either 
at our accuracy. Thus the result 1(75(1 from the GSM regime remains valid in the Aoki regime, in agreement with 
Refs. BE]!. 

The other quantities which are substantially changed in the Aoki regime are the parity violating matrix elements. 
These are 0(a) in the GSM regime, but become of 0(1) in the Aoki regime. This is indicated by the factors of a 2 /M' 
in the W' contributions in Eqs. 1(91196(1 . In fact, it is straightforward to see that we control these quantities only at 
O(l), and not at O(a). This simplifies the calculation since only the pole terms give O(l) contributions. The results 
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are (with a = 1, 2): 

M 1 sin e 

<7r a (p 2 )|A°|7r 3 (pi)> = g M 2 + 0(a) , (109) 

(7r a (p2)|i3|7r a (pi)) = J<7r 3 (p2)|iJ|7r3(pi)> = + °( a ) > ( n0 ) 

<^(p 2 )|P>3bi)) = -»i?o f! Sm 2 £ + 0(a) , (111) 

MP2)\P 3 K(Pi)) = UirM\P 3 \irM) = -iB f*^ + 0(a) . (112) 
3 q z + m£ 3 

We have used Eq. I|102|) to simplify the expressions. Note that the choice of pion mass in the denominators makes an 
O(l) difference to these quantities. 

The expressions in this section all simplify at maximal twist. As already noted, we would then have e = in 
addition to c = and s = 1. Thus the O(a) corrections to the physical quantities in H104I107[) vanish, so that all are 
given simply by their LO continuum xPT values. The charged pion masses are similarly given by the LO continuum 
result, = M'. The pion mass-squared splitting does differ from the continuum, but this difference is of 0(a 2 ). 
Finally, the parity violating matrix elements become of 0(a) again, as in the GSM regime, although we cannot predict 
them at the order we are working. 



We now return to the issue raised in Ref. |2l|. The authors argue that automatic O(a) improvement at maximal 
twist does extend into the Aoki regime, but only if one uses a twist angle based on a critical mass which includes the 
O(a) offset, to' — m. They point out that one of the definitions of twist angle that accomplishes this is what we call 
uj as defined by Eq. (and which they call lowt)- We fully agree on the need to use such a definition, not only in 
the Aoki regime, but also in the GSM regime, as we discussed in Sec. [HI 

As for automatic 0(a) improvement, our results support the proposal of Ref. |2l|, although with one caveat. We 
argued that, in the Aoki regime, the condensate did not vary smoothly as the masses are varied. This is true in 
general, but not at maximal twist, where the condensate remains fixed, with uo m = tt/2 within an uncertainty of 0(a). 
Indeed, as we have already noted, at maximal twist physical quantities are described in the Aoki phase by continuum 
%PT up to 0(a 2 ). If we start in the Aoki regime at u> = n/2, and remain on the [i axis as a is reduced (so that we 
move into the GSM regime and ultimately to the continuum limit) we expect that physical quantities, at fixed quark 
masses, will extrapolate smoothly to the continuum with errors quadratic in a. 

The caveat is that the argument will fail if, as [i is reduced, one encounters a phase boundary. This is expected 
to happen in one of the two possible scenarios for the phase diagram [l3L ITiL [l5 | . In this scenario, in which W' > 0, 
the end-point of the phase boundary occurs at fi = 32d 2 W^'// 2 , where m^ 3 = 0. Values of fi greater than this are 
allowed, but for smaller values the condensate changes rapidly (by an amount of O(l) while the quark mass changes 
by 0(a 2 )), and continuum xPT expressions fail. In the other scenario (W < 0) we see no need to impose a lower 
limit on /i — unless \W'\ is unnaturally small, yet higher order terms in the chiral potential will not change the phase 
structure. 

It is worth noting the accuracy with which the critical mass must be determined to stay at lu — ir/2 within an error 
of O(a). This is the accuracy that is required to automatically remove the O(a) terms. In the GSM regime, where 
/i ~ aAq CD , this requires knowing to' to an absolute accuracy of 0(a 2 ). In the Aoki regime, however, the required 
absolute accuracy decreases to 0(a 3 ). This sounds like a difficult goal, but, in fact, as we discussed in the summary 
of the previous section, it is attainable even by doing simulations in the GSM regime. 



V. CONCLUSION 



In this paper we have used effective field theory methods to study the discretization errors in tmLQCD. We have 
presented results for a number of pionic quantities that can be calculated in lattice simulations, and studied different 
possible definitions of the twist angle. Perhaps the most interesting quantities to calculate are the difference to — ujp 
(which indicates the size of the O(a) uncertainty in the twist angle), the splitting m 2 3 — to 2 ^ (which indicates the size 
of 0(a 2 ) flavor breaking, but parity conserving, quantities) and the axial and pseudoscalar form factors of the pion 
(which are the simplest examples of unphysical quantities). To our knowledge, none of these quantities have been 
calculated to date, and we urge that they be considered in the future. 

In the introduction, we raised a number of issues concerning tmLQCD; now we can comment on what we have learnt 
about these. Concerning the impact of flavor and parity breaking, we have seen how this contributes to unphysical 
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quantities at 0(a), and in one case, to physical quantities at 0(a 2 ). This case is the pion mass splitting, where we 
confirm the interesting result of Ref. [l4( that knowledge of the sign of the splitting allows one to predict the nature 
of the phase structure in the Aoki region. Nevertheless, at NLO in tm^PT, the pions in the loops are degenerate, 
and there is no flavor breaking in decay constants or vector and scalar form factors. For staggered fermions, where 
the "taste" breaking between pions is also of 0(a 2 ), it has been found essential to include this breaking in loop 
contributions in order to make good fits to the chiral behavior j^lj. If this turns out to be true also for tmLQCD 
(which depends on the size of the as yet unmeasured flavor breaking), then the %PT calculations will need to be 
extended beyond NLO. One possibility which we are investigating is to use a power counting like that in Ref. [3l]], so 
that a full NNLO calculation is not needed. 

We also asked which quantities are good indicators of the size of discretization errors. One answer is that it may 
be easier to use the flavor-parity violating quantities (i.e. the difference between different definitions of twist angle, 
and the axial and pseudoscalar form factors), since these are of 0(a) (even at maximal twist) and thus easier to 
calculate. In more general terms, studying these quantities in simulations would allow a more thorough test of our 
understanding of tmLQCD. 

Another interesting question is how we can test the reliability of tm%PT at the order we are working. Our results, in 
fact, contain a number of predictions that are valid at NLO but not at higher order. In particular, as we have outlined 
in the summary subsection of Sec. lIIII most of the flavor-parity violating form factors can be predicted once the other 
quantities that we have discussed have been calculated. Such tests are important for establishing the credibility of 
the chiral and continuum extrapolations that must ultimately be done. 

We have discussed extensively the issue of the smallest quark mass that can be used without invalidating automatic 
0(a) improvement at maximal twist. We have argued that one can certainly work in the GSM regime (m q ~ qAq CD ) 

and also at least part-way into the Aoki regime (m q ~ a 2 Aq CD ). This is only true, however, if one uses the appropriate 
definition of twist angle. Allowable choices are that determined from the vanishing of the coupling of the bare axial 
current in the twisted basis to the charged pion, or from the vanishing of the coupling of the pseudoscalar density 
to the neutral pion. Alternatively one could fix the twisted mass /i and vary the untwisted mass until the pion 
masses (either charged or neutral) reach their minimum values Jllj . What is potentially problematic is to determine 
the critical quark mass by extrapolating the squared pion mass to zero along the untwisted Wilson axis using only 
moderately small quark masses. In the GSM regime this may not, in practice, give the critical mass with the accuracy 
needed (errors of absolute size 0(a 2 )) to ensure 0(a) improvement at maximal twist. 6 Furthermore, this method does 
not apply in the Aoki regime, where it gives a determination of the critical mass with an absolute error of 0(a 2 ), 
whereas the required accuracy is 0(a 3 ). 

One spin-off from our calculation is a method for determining the new LECs that enter when one incorporates 
discretization errors. This both allows a direct measure of their size, and can be used a posteriori to correct results 
obtained with untwisted Wilson fermions. Furthermore, our calculations extend the results that are available from 
XPT applied to untwisted Wilson quarks 0,^3 t° several new quantities. 

Many of the quantities we have considered will be difficult to calculate in lattice simulations because they involve 
quark-disconnected contractions. An important practical question to consider is whether, by using partially quenched 
tmLQCD and its corresponding xPT J32j, one can make predictions for the quark-connected and disconnected con- 
tributions separately. It will also be interesting to extend the calculations to the case of a non-degenerate doublet of 
quarks [33|. 
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This point may explain the "bending phenomenon" observed at masses below m q ~ &Aq CD in Ref. 1 1(1 . Such bending is expected in 
most quantities if one approaches the untwisted Wilson axis along a line parallel to the twisted mass axis. It would be interesting to fit 
the results from this work under this hypothesis, using the formulae we have provided. 
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